Gated Diffusion-controlled Reactions

The binding and active sites of proteins are often dynamically occluded by motion of the nearby polypeptide. A variety of theoretical and computational methods have been developed to predict rates of ligand binding and reactivity in such cases. Two general approaches exist, "protein centric" approaches that explicitly treat only the protein target, and more detailed dynamical simulation approaches in which target and ligand are both treated explicitly. This mini-review describes recent work in this area and some of the biological implications

Molecular Basis of Ligand Dissociation in β-Adrenergic Receptors

The important and diverse biological functions of β-adrenergic receptors (βARs) have promoted the search for compounds to stimulate or inhibit their activity. In this regard, unraveling the molecular basis of ligand binding/unbinding events is essential to understand the pharmacological properties of these G protein-coupled receptors. In this study, we use the steered molecular dynamics simulation method to describe, in atomic detail, the unbinding process of two inverse agonists, which have been recently co-crystallized with β1 and β2ARs subtypes, along four different channels. Our results indicate that this type of compounds likely accesses the orthosteric binding site of βARs from the extracellular water environment. Importantly, reconstruction of forces and energies from the simulations of the dissociation process suggests, for the first time, the presence of secondary binding sites located in the extracellular loops 2 and 3 and transmembrane helix 7, where ligands are transiently retained by electrostatic and Van der Waals interactions. Comparison of the residues that form these new transient allosteric binding sites in both βARs subtypes reveals the importance of non-conserved electrostatic interactions as well as conserved aromatic contacts in the early steps of the binding process

Effects of Macromolecular Crowding on Protein Conformational Changes

Many protein functions can be directly linked to conformational changes. Inside cells, the equilibria and transition rates between different conformations may be affected by macromolecular crowding. We have recently developed a new approach for modeling crowding effects, which enables an atomistic representation of “test” proteins. Here this approach is applied to study how crowding affects the equilibria and transition rates between open and closed conformations of seven proteins: yeast protein disulfide isomerase (yPDI), adenylate kinase (AdK), orotidine phosphate decarboxylase (ODCase), Trp repressor (TrpR), hemoglobin, DNA β-glucosyltransferase, and Ap4A hydrolase. For each protein, molecular dynamics simulations of the open and closed states are separately run. Representative open and closed conformations are then used to calculate the crowding-induced changes in chemical potential for the two states. The difference in chemical-potential change between the two states finally predicts the effects of crowding on the population ratio of the two states. Crowding is found to reduce the open population to various extents. In the presence of crowders with a 15 Å radius and occupying 35% of volume, the open-to-closed population ratios of yPDI, AdK, ODCase and TrpR are reduced by 79%, 78%, 62% and 55%, respectively. The reductions for the remaining three proteins are 20–44%. As expected, the four proteins experiencing the stronger crowding effects are those with larger conformational changes between open and closed states (e.g., as measured by the change in radius of gyration). Larger proteins also tend to experience stronger crowding effects than smaller ones [e.g., comparing yPDI (480 residues) and TrpR (98 residues)]. The potentials of mean force along the open-closed reaction coordinate of apo and ligand-bound ODCase are altered by crowding, suggesting that transition rates are also affected. These quantitative results and qualitative trends will serve as valuable guides for expected crowding effects on protein conformation changes inside cells

Direct observation of topoisomerase IA gate dynamics

Type IA topoisomerases cleave single-stranded DNA and relieve negative supercoils in discrete steps corresponding to the passage of the intact DNA strand through the cleaved strand. Although type IA topoisomerases are assumed to accomplish this strand passage via a protein-mediated DNA gate, opening of this gate has never been observed. We developed a single-molecule assay to directly measure gate opening of the Escherichia coli type IA topoisomerases I and III. We found that after cleavage of single-stranded DNA, the protein gate opens by as much as 6.6 nm and can close against forces in excess of 16 pN. Key differences in the cleavage, ligation, and gate dynamics of these two enzymes provide insights into their different cellular functions. The single-molecule results are broadly consistent with conformational changes obtained from molecular dynamics simulations. These results allowed us to develop a mechanistic model of interactions between type IA topoisomerases and single-stranded DNA

The Effect of Macromolecular Crowding, Ionic Strength and Calcium Binding on Calmodulin Dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM's target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.Comment: Accepted to PLoS Comp Biol, 201

Peptide Folding in Cellular Environments: A Monte Carlo and Markov Modeling Approach

Steric interactions with surrounding macromolecules tend to favor the compact native state of a globular protein over its unfolded state. However, in experiments conducted in cells and concentrated protein solutions, both stabilization and destabilization of proteins have been observed, compared to dilute-solution conditions. Therefore, in order to understand the effects of surrounding macromolecules on protein properties such as stability, there is a need for computational modeling beyond the level of hard-sphere crowders. Here, we discuss some recent exploratory studies of peptide folding in the presence of explicit protein crowders, carried out by us using an all-atom Monte Carlo-based approach along with an implicit solvent force field. For interpreting the simulation data, time-lagged independent component analysis and Markov state modeling are used